A new route to Explosive Percolation 
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The biased link occupation rule in the Achlioptas process (AP) discourages the large clusters to 
grow much ahead of others and encourages faster growth of clusters which lag behind. In this paper 
we propose a model where this tendency is sharply reflected in the Gamma distribution of the cluster 
sizes, unlike the power law distribution in AP. In this model single edges between pairs of clusters of 
sizes Si and Sj are occupied with a probability oc (siSj) a . The parameter a is continuously tunable 
over the entire real axis. Numerical studies indicate that for a < a c the transition is first order, 
q c = for square lattice and a c = —1/2 for random graphs. In the limits of a = — oo,+oo this 
model coincides with models well established in the literature. 
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A first order transition in critical phenomena is char- 
acterized by an abrupt jump in the order parameter 
An infinitesimal increase of the control variable leads to 
a spurt of activity on the global scale as in a catas- 
trophic process. The percolation transition describes a 
geometrical phase transition between ordered and disor- 
dered phases of random resistor network, binary alloys, 
forest fires, galaxies etc. which are continuous transi- 
tions [2-4]. On the other hand first order transition is 
observed in Bootstrap Percolation [f| modeling the com- 
petition between exchange and crystal field interaction 
in some magnetic materials. Here, depending on the lat- 
tice structure, it may happen that culling of even a single 
spin evacuates a globally connected system in a recursive 
process. 

Recently it is observed that a biased link occupation 
rule, known as the Achlioptas process |7|, on the Erdos- 
Renyi random graph [|| leads to a first order transition. 
Here a pair of edges (ij) and (hi) are randomly selected 
between clusters of sizes Sj, Sj and Sk, sz respectively. 
The edge with the smaller value of the products SiSj and 
Sk-si is occupied; when equal, one edge is selected ran- 
domly. Repeated application of this rule for every edge 
delays growth of the largest cluster. Therefore the biased 
edge occupation rule in AP discourages a single cluster 
to grow much ahead of other clusters and at the same 
time it encourages faster growth of clusters which lag 
behind. One therefore expects that at a certain stage 
the cluster size distribution should have a characteristic 
scale. Finally nearly equal size clusters are linked to- 
gether by additional few edges (of zero measure) resulting 
an abrupt global connection as in first order transition, 
hence the name 'Explosive Percolation' (EP) [7], trig- 
gering an explosion of scientific activity. Later, details 
of AP on square lattice Q and random graph § have 
been studied. Also AP on scale- free networks observes 
a transition from a continuous to a first order transition 
Hamiltonian formulation of AP is studied in 
13j . first order transition has been observed in a cluster 
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FIG. 1: A small graph of 5 nodes and 2 links. The next link 
is placed using the probabilty 7Ty as defined in the text. 



aggregation process L4| and also in human protein ho- 
mology network 15] . However, very recently, Costa et. 
al. has claimed that the nature of transition in Explosive 
Percolation is indeed continuous |16j . Opposite bias can 
also be given by occupying the edge with larger product 
and this percolation transition remains continuous. 

In this paper we argue that it is possible to attain EP in 
a generic fashion by introducing a continuously tunable 
parameter and selecting only a single edge at a time. At 
each step in the process we consider all currently vacant 
edges (if we wish to avoid edges internal to existing clus- 
ters, the loop-less condition defined below, we consider 
only those edges which connect clusters). Let the sizes of 
the respective clusters to be connected by the edges be 
denoted Sj and Sj . Each of these vacant edges is given a 
weight (siSj) a and the sum over all the weights is calcu- 
lated as a normalization constant W . Then one of these 
vacant edges is chosen with probability (siSj) a /W and 
occupied. This procedure is then repeated for the next 
edge addition. 

For example, let us consider a N = 5 nodes graph 
(Fig. 1). Let at an intermediate stage there be only two 
links in this graph, connecting nodes 1 and 5 and nodes 
3 and 4. Therefore the graph has three components of 
sizes 2, 2, 1 and has 8 vacant edges: (1,2), (1,3), (1,4), 
(2,3), (2,4), (2,5), (3,5) and (4,5). Now, to place the 
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FIG. 2: Percolation with a — 5* — oo on square lattice (Z = 24) and with loop-less condition, (a) Only links are deposited (b) 
clusters of size up to 4 (c) clusters of size up to 8 and the (d) the spanning tree. Similarly percolation with a — > +oo are: (e) 
growth starts from a single link (f) cluster of N/4 bonds (g) cluster of N/2 bonds (h) the spanning tree cluster. 



third link in this graph we calculate weights associated 
with every remaining vacant edge of the graph. For some 
specific pre-assigned value of a the associated weights 
are: w 12 = (2.1)", w 13 = (2.2) Q , w u = (2.2) Q , w 23 = 
(1.2)°*, w 2i = (1.2)", w 25 = (1.2)« w 35 = (2.2)", w i5 
= (2.2)". The sum of all these weights: W = Sw 13 
= wi2+wi3+wi4+W23+W24+W25+if35+it;45 and the edge 
occupation probablities are 7Ty = uiij/W. One vacant 
edge is then randomly selected with probability Try and 
is occupied. This completes the procedure to occupy one 
edge which is then repreated. 

Here a being a continuously tunable parameter vary- 
ing over the entire real axis from — oo to +oo. This is 
the crucial difference between AP and our model. In 
AP both links are picked up with uniform probability as 
in random graph, i.e., each link connects clusters with 
probability oc SiSj. In comparison, in our case links are 
picked up using cluster size dependent probabilities. For 
a < 0, edges with smaller products are preferentially se- 
lected which delays the transition and in the following we 
present numerical evidence that for all a < a c the per- 
colation transition is first order. For a > edges with 
larger products are occupied with higher probabilities. 
Our model is completely distinct from AP, i.e., AP can- 
not be retrieved for any value of a except in two limits. 
AP can be extended by selecting the minimal product of 
n randomly selected edges. In the limit of n — > oo the 
edge which has the globally minimal (maximal) value of 
the product is only selected. Only these two limits iof ex- 
tended AP correspond to a — > — oo(+oo) of our model. 

The 'sequential time' t is defined as the number of 
links occupied and the link densities are r = t/N for 
random graph and p = t/(2N) for square lattice with 
periodic boundary condition. We consider percolation 
with (without) loops where links are allowed (forbidden) 
to connect two nodes of the same cluster. The ordinary 



bond percolation corresponds to a = 0. When 'loop-less' 
condition is imposed this is called 'loop-less percolation' 
fl7j ] with p c = (7 — 3-\/3)/4 on square lattice 0. Since 
random graph is infinite dimensional loop-less condition 
has no effect and r c = 1/2. 

Next, the a — > — oo loop-less case is considered. Here 
at any stage one edge is randomly occupied only from 
the subset of all vacant edges for which the SiSj val- 
ues are minimum. In the beginning all N nodes are 
isolated and the minimal product (MP) is 1. The MP 
condition is essentially a hard core repulsion and there- 
fore successive links are placed without overlapping with 
the previous ones. This process is known as the "ran- 
dom sequential adsorption" (la ]. Gradually the system 
reaches a jamming limit when no more link can be placed 
(Fig. 2(a)). The jammed state of random graph has N/2 
links but there are isolated nodes and links in square lat- 
tice. At this stage the MP jumps to 4 for random graph 
but to 2 for square lattice. Consequently the next se- 
quence of links connect pairs of links in random graph 
but single sites with links in square lattice. Eventually 
the vacant edges with this MP are again exhausted at a 
second jammed state when the MP is further enhanced 
and this process continues in a series of jammed states 
with discrete jumps of MPs (Fig. 2(b-c)). The sequence 
of MPs for square lattice are 1, 2, 4, 6 etc. with jam- 
ming densities 0.227(1), 0.273(1), 0.351(1), 0.370(1) re- 
spectively which approach the p c = 1/2, a spanning tree 
configuration (Fig. 2(d)). In contrast the jamming den- 
sities in random graph are (s — l)/s for clusters of sizes 
s = 2,4,8,16,... etc. giving r c = 1. The largest cluster 
size when plotted with the link density gives a step func- 
tion at r = 1 in iV - >oo,a perfect first order transition 
also observed in [9(. Similarly for a — > +oo vacant edges 
with maximal values of SiSj are only occupied (Fig. 2(e- 
h)). Once the first link is placed, the subsequent links 
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FIG. 3: (Color online) (a) Order parameter C(p,a) with link 
density p for a = 1/2, 0.2, 0.1, 0, -0.1, -1/4, -1/2, -1, -2, -3, 
-4 and -5, a values decreasing from left to right, for square 
lattice of L — 512. (b) The asymptotic values of the gap 
5(a), the percolation threshold p c (a) — 1/2 and the largest 
jump g(a) of the order parameter plotted with a. 



get connections to this cluster only, therefore only one 
cluster grows. In random graph a new node is connected 
to a randomly selected node of the growing cluster. This 
is the 'model A' network of [l9| with an exponentially de- 
caying degree distribution. On square lattice the growth 
is limited to only the surface bonds leading to the well 
known 'Bond Eden Tree' [13] ■ Without loop- less restric- 
tion all vacant edges within the cluster are first occupied 
and bulk of the cluster grows like a perfect crystal. 

Now we study the percolation process with loops for 
finite value of a. Fig. 3(a) plots the order parameters 
C[p,a) = s m (p,a)/N, s m being the size of the largest 
cluster, for square lattice with size L = 512; N = L 2 . 
As a decreases from zero the order parameter increases 
more rapidly and the whole curve shifts to the larger p 
values indicating increasingly delayed transitions. For 
—5 < a < —1 the curves are nearly vertical, C(p, a) 
jumps from to 1 within a short interval of p. However 
within — 1 < a < 0, the curves are relatively smooth but 
gradually become steeper as N increases. In the following 
we present evidence that for all a < a c — the perco- 
lation transitions are likely to be first order transition, 
with p c > 1/2. 

To determine the order of transition and estimate the 
percolation threshold we define: to is the latest time with 
s m < N 1 / 2 and t\ is the earliest time s m > kN; k = 1/2 
(random graph) [tJ and 0.9 (square lattice). A gap A = 
t\ — tg in N — > oo limit can distinguish between two kinds 
of transitions: 5(a) = limjv->oo ^°^ N ^ = constant > 
for a continuous transition and for a first order transi- 
tion. Average values of A (a, N) /N are then extrapolated 
to N — !• oo limit as N~@( a \ Different /3(a) values for 
an a are tried and the one with the minimal Standard 



FIG. 4: (Color online) (a) For AP on random graph, the 
maximal jump g(N) in the order parameter has been plotted 
on a double logarithmic scale with graph size N. The slope is 
-0.06475. (b) The maximal jump g(N) is then plotted against 
jy-0-06475 Qn & jj near sca i e ^he continuous line is a straight 
line fit of this data which is then extrapolated to N — >• oo to 
meet g(N) axis at -0.000516. 



Error 21j is selected. The 5(a) values are nearly zero 
for —5 < a < —1/4 indicating the first order transition 
(Fig. 3(b)). There is a gradual increase of 5(a) within 
— 1/4 < a < due to the finite size of systems simulated. 
However 5(0) w 0.101 is obtained for ordinary percola- 
tion. 

We also studied in detail the extent of jump in the or- 
der at the transition point. For that we calculated the 
incremental changes in C due to the occupation of an ad- 
ditional link and keep track of these changes as links are 
occupied one by one. The average value of the largest 
jump AC m due to addition of only a single link and 
the corresponding link density r m have been calculated. 
In the thermodynamic limit of N — > oo the asymptotic 
jump length is: g = linijv->.oo g(N) = liniAr^oo AC m (N). 
We expect g > for a discontinuous transition and 
g = for a continuous transition. We first calculate this 
quantity for the AP on random graph. This data has 
been shown in Fig. 4 and g(N) is observed to vary as: 
g(N) - iv -0 064771 which implies that g(N) approaches 
to g = as N — > oo. We believe this result is consis- 



tent with the claim in |16| that the AP is a continuous 
transition. 

For our problem the asymptotic jump g(a) = 
lim^v-i-oo AC m (a, N) for different a are plotted in Fig. 
3(b). For ordinary bond percolation g(a — 0,N) — > 
as iV"^/ 2 - 1 , with d f = 91/48, the fractal dimension of 
the incipient infinite percolation cluster Q . However for 
a < the g(a) jumps to 0.16 at a = —0.05 and then 
gradually increases to 1/2 as a -> -oo. The percolation 
thresholds p c (N) are estimated by the average values of 
r m and t\ giving approximately same asymptotic value 
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FIG. 5: (Color online) (a) Order parameter C(r,a) with link 
density r for a = 1/2, 1/4, 0, -0.1, -1/4, -1/2, -1, -2, -3, -4 and 
-5, a values decreasing from left to right, for random graph 
of N = 4096. (b) The asymptotic values of the gap 8(a), the 
percolation threshold r c (a) — 1/2 and the largest jump g(a) 
of the order parameter plotted with a. 

for p c (a). The p c (0) = lim;v-voo p c (0, N) ->• 1/2 limit is 
approached as N~ x / 2v , = 4/3 [3j being the correlation 
length exponent for the ordinary percolation. On the 
other hand for a < 0, p c (c*0 increased continuously with 
a (Fig. 3(b)) and tends to unity as a — > — oo. 

In Fig. 5(a) we show C(r, a) vs. r for random graph. 
Last five curves for a < — 1 almost coincide with the 
vertical line at r — 1 implying first order transition with 
r c = 1 in this range. Curves are smoother for — 1 < a < 
0. In Fig. 5(b) we show r c (a) tends to 1/2 rapidly as 
a 0. The jump g(a) is almost zero for— 1/4 < a < 
but then it rapidly increases and tends to 1/2 as a -> 
— oo. The gap 8(a) is almost zero for a < —1/2 and 
then slowly increases to ps 0.193 for random graph. This 
result indicates that for random graph a c = —1/2. 

The approach to first order transition in our system 
is completely different from AP. Comparison is made 
in the tendency of supressing growth of large clusters 
and enhancing growth of small clusters as reflected in 
following three quantities. In Fig. 6(a) we show the 
scaled cluster size distribution P(s,p) for square lattice 
at a = —1. The individual distributions has strong de- 
pendence on Ap — p c — p but not significantly on L. 
The best scaling form is P(s,p)Ap~ v ~ Q(sAp^) with 
r] = ( = 1.77(5). The scaling function fits to the Gamma 
distribution Q(x) ~ x a exp(-bx) with a = 0.93(10) and 
b = 4.09(10). P(s,p) for a < -1 also fits to Gamma 
distributions with different values of a and b. In addition 
we have checked the finite size scaling analysis as well. 
Right at the percolation point the cluster size distribu- 
tion P(s, P c , L)L 9 scales excellent with s/L 09 and the 
scaling function fits very well to a Gamma distribution. 
Here the percolation point is determined by the maximal 
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FIG. 6: (Color online) (a) Finite-size scaling of the cluster 
size distribution P(s,p) of square lattice at a = — 1 (with 
p c = 0.549(2)) for L= 64 and 128; p = 0.48, 0.50 and 
0.52. The solid line is a fit to the Gamma distribution, (b) 
The binned average cluster size distribution P(s,p) of AP on 
square lattice for L = 1024; p = 0.49, 0.51, 0.52 and 0.525 
(from left to right). 



jump in the order parameter. 

In comparison P(s,r) in AP follows a power law as 
in continuous transition and the exponent approaches to 
R3 —2.11(2) as r — > r c . As r c is approached, a hump 
appears in the tail implying an enhanced population of 
large clusters. At percolation point these clusters get 
linked giving rise to a sudden jump of the size of the 
largest cluster. We believe this is due to random selec- 
tion of edges in AP the power law distribution masks 
the Gamma distribution; the exponent is raised to -2.11 
from -5/2 in random graph but retains -2.06(1) for square 
lattice (Fig. 6(b)) compared to 187/91 for ordinary per- 
colation in 2d also observed in (l2j . 

Since in first order transition there are many clusters 
whose sizes are nearly same we calculate the average ra- 
tio lZ(p) — (s m / ' s nm ) where s nm is the size of the sec- 
ond largest cluster. For first order transition we expect 
7Z(p) to maintain a uniform variation with p where as it 
should increase smoothly beyond the percolation thresh- 
old for a continuous transition. Fig. 7(a) shows TZ(p) for 
square lattice with L = 128 and AP on square lattice. 
For all a < the plots are indeed horizontal lines up to 
p = Pc(a>). It implies that s m and s nm have comparable 
values till the system reaches p c and the lines stop shortly 
after p c since at percolation threshold the largest cluster 
covers the entire system. In comparison lZ(p) for a = 
and AP on square lattice grows continuously to very large 
values indicating that the largest cluster grows by swal- 
lowing small clusters and it continues far beyond p c . A 
third quantity is the decay of the fraction ni(p, a) of iso- 
lated nodes in the system. For the first order transition, 
the ni(p, a) decays very fast and vanishes at or before p c 
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FIG. 7: (Color online) (a) Ratio lZ(p, a) of the sizes of 
largest cluster and the next largest with p/p c (a) for a — 
—3,-2,-1,-1/2,0 and AP on square lattice from bottom 
to top. For more visibility y axes have been multiplied by 1, 
2, 4, 8, 16 and 32 respectively, (b) Fraction n\(jp,ct) of iso- 
lated nodes with p/p c (a) for the same a values (from right to 
left) as in (a). 



for all a < (Fig. 7(b)). It is also observed that as a 
decreases ni(p, a) vanishes even faster, e.g., as a — > — oo, 
ni(p, — oo) decreases linearly to zero at p = 0.273. In 
comparison n\(p) decays as exp(— 2r) in random graph 



and exp(- 



.1.17 



) in the AP, though on square lattice they 



behave different from our model (Fig. 7(b)). We ar- 
gue that the cluster size distribution P(s,p), ratio TZ(p) 
of two largest clusters and the density of isolated nodes 
ni (p) reflect distinguishing signatures of our model with 
AP at, beyond and before the transition. 

Finally we like to mention that numerical study of our 
model is CPU intensive since for a ^ the CPU ~ A^ 3 for 
random graph and ~ A^ 2 for square lattice in comparison 
to ~ A" in AP. Therefore we could study up to L = 512 for 
square lattice and A" = 4096 for random graph, though 
the extent of our computational efforts are comparable 
to other studies in the literature. 

To summarize, we have studied a new route to Explo- 
sive Percolation where single edges between pairs of clus- 
ters are randomly occupied with probabilities oc (siSj) a . 
From our numerical studies we present evidence that for 
all a < a c the transitions are discontinuous first order 
transitions. We obtain a c = for the square lattice but 
-1/2 for random graphs. Our model is completely differ- 
ent from AP apart from the limits of a — > — oo and +oo 
where it coincides with similar limits of the generalized 



AP. The effect of the biased occupation rule discourag- 
ing growth of large clusters and encouraging growth of 
small clusters is distinctly visible in the Gamma distri- 
bution of the cluster sizes in our model. This is in con- 
trast to power law distribution in AP like the continuous 
transition which we believe is due to random selection 
of edges. In addition we define two quantities 7Z(p) and 
n-i(p) which distinguishes our model from AP. 
manna@bose.res.in 



[2] 
[3] 

[4] 

[o] 

[6] 

[7] 

[8] 
[9] 

[10] 

[11] 
[12] 
[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 
[20] 



[1] H. E. Stanley, Introduction To Phase Transitions And 

Critical Phenomena, Oxford University Press, 1971. 
2] G. Grimmett, Percolation, Springer, 1999. 
3] D. Stauffer and A. Aharony, Introduction to Percolation 

Theory (Taylor & Francis, London, 1994). 
4] P. W. Kasteleyn and C. M. Fortuin, J. Phys. Soc. Japan 
26, 11 (1989). 

5] J. Chalupa, P. L. Leath and G. R. Reich, J. Phys. C, 

12, L31 (1979). 
6] P. Erdos and A. Renyi, Publ. Math. Debrecen, 6, 290 
(1959). 

7] D. Achlioptas, R. M. D'Souza, and J. Spencer, Science 

323, 1453 (2009). 
8] R. M. Ziff, Phys. Rev. Lett. 103, 045701 (2009). 
9] E. J. Friedman A. S. Landsberg, Phys. Rev. Lett. 103, 
255701 (2009). 

0] Y.S. Cho, J. S. Kim, J. Park, B. Kahng and D. Kim, 

Phys. Rev. Lett. 103, 135702 (2009). 
1] F. Radicchi and S. Fortunato, Phys. Rev. Lett. 103, 
168701 (2009). 

2] F. Radicchi and S. Fortunato, Phys. Rev. E 81, 036110 
(2010). 

3] A. A.Moreira, E. A. Oliveira, S. D. S. Reis, H. J. Her- 
rmann and J. S. Andrade Jr. Phys. Rev. E 81, 040101 
(2010). 

4] Y. S. Cho, B. Kahng and D. Kim, Phys. Rev. E 81, 
030103 (2010). 

5] H. D. Rozenfel d, L. K. Gallos and H. A. Makse, 

larXiv:0911.4082l 
6] R. A. da C osta, S. N. Dorog ovtsev, A. V. Goltsev, J. F. 

F. Mendes. larXiv: 1009.25341 
7] S. S. Manna and B. Subramanian, Phys. Rev. Lett. 76, 
3460 (1996). 

3] S. S. Manna and N. M. Svrakic, J. Phys. A., 24, L671 
(1991). 

9] A.-L. Barabasi and R. Albert, Science, 286, 509 (1999). 
0] D. Dhar and R. Ramaswamy, Phys. Rev. Lett. 54, 1346 
(1985), S. S. Manna and D. Dhar, Phys. Rev. E. 54, 
R3063 (1996). 

[21] mathworld . wolfram . com/LeastSquaresFitting . html. 



